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Abstract 

Many special functions are solutions of first order linear systems y'n{x) = 
an{x)yn{x) + dn{x)yn-i{x), y'n-i{x) = 6„(x)y„_i(x) + en{x)yn{x). We ob- 
tain bounds for the ratios and the logarithmic derivatives 
of yn{x) for solutions of monotonic systems satisfying certain initial condi- 
tions. For the case dn{x)en{x) > 0, sequences of upper and lower bounds 
can be obtained by iterating the recurrence relation; for minimal solutions 
of the recurrence these are convergent sequences. The bounds are related to 
the Liouville-Green approximation for the associated second order ODEs as 
well as to the asymptotic behavior of the associated three-term recurrence 
relation as n — t- +oo; the bounds are sharp both as a function of n and 
X. Many special functions are amenable to this analysis, and we give sev- 
eral examples of application: modified Bessel functions, parabolic cylinder 
functions, Legendre functions of imaginary variable and Laguerre functions. 
New Turan-type inequalities are established from the function ratio bounds. 
Bounds for monotonic systems with dn{x)en{x) < are also given, in par- 
ticular for Hermite and Laguerre polynomials of real positive variable; in 
that case the bounds can be used for bounding the monotonic region (and 
then the extreme zeros). 
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1 Introduction 

Many special functions, and in particular functions of hypergeometric type, 
satisfy first order differential systems of the form 



For the particular case of modified Bessel functions sharp bounds for 
function ratios and logarithmic derivatives y'^{x) / yn{x) , as 

well as Turan-type inequalities were recently obtained in [T^; the key ingre- 
dient in the analysis was the study of the qualitative behavior of the solutions 
of the Riccati equation satisfied by hn{x) = together with 

the application of the three-term recurrence relation. 

Ratios of Bessel functions appear in a great number of applications, 
particularly as parameters of certain probability distributions (see, for in- 
stance, the examples mentioned in |1]). Other special function ratios are 
important in applications. In particular, parabolic cylinder ratios appear in 
the study of Ornstein-Uhlenbeck processes (see, for instance [2]), and other 
special function ratios (Whittaker, Legendre, Gauss hypergeometric func- 
tions) play similar roles as well OHUS]. In all these applications, a common 
characteristic is that the functions are real and the variables lie inside a 
monotonic region (region free of zeros). 

In this paper, we put the ideas of [T] in a more general context and we 
analyze the qualitative behavior of the Riccati equation associated to the 
ratio hn{x) = ynix)/y„^i{x), 



y'^{x) = an{x)yn{x) + (i„(x)y„_i(x), 

y'n^li^) = hn{x)yn-l{x) + en{x)yn{x). 
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in the general case in which the quadratic equation 

enix)Xnixf + {bn{x) - a„(x))An(x) - = (3) 

has two distinct real roots A^(x). This case corresponds to monotonic sys- 
tems, with solutions which have one zero at most. As we will see, if the 
functions A^(x) are monotonic, they are bounds for the ratios hn{x) satis- 
fying certain initial value conditions. 

We will discuss in detail bounds for systems with dn{x)en{x) > 0. Up- 
per and lower bounds become available which are accurate approximations 
for large parameters and/or variable and with bounded error (upper and 
lower bounds are available). The bounds given are related to the Liouville- 
Green approximation for the associated second order ODE as well as to 
the asymptotic behavior of the associated three-term recurrence relation. 
This explains why the bounds become sharper as the variable x and/or the 
parameter become large. 

We provide several examples (section 14. ip : modified Bessel functions, 
parabolic cylinder functions. Associated Legendre functions of imaginary 
argument and Laguerre functions of negative argument. From these bounds, 
a good number of new Turan-type inequalities are obtained. Turan-type 
properties for special functions have received a considerable attention in 
recent years; just to cite five different groups of researchers, we mention 
[3 [71 |8l [9l d] (see also references cited therein) . 

We also give two examples of applications of the methods for the case 
dn{x)en{x) < (section l4.2p and use these results for bounding function 
ratios for Laguerre and Hermite polynomials in the real axis (but outside the 
oscillatory region). These bounds can be used for bounding the oscillatory 
region and, therefore, for bounding the extreme zeros. 

In addition to direct applications in several areas, particularly in statis- 
tics and stochastic processes, the bounds on function ratios have implica- 
tions in the construction of numerical algorithms. These techniques provide 
bounds for the region of computable parameters of a given function within 
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the overflow and underflow limitations, and they also provide bounds for 
the condition numbers of the functions (see section 14.1.21 for the case of 
Parabolic Cylinder Functions). Additionally, as discussed for the particular 
case of modified Bessel functions [l], the bounds are useful for accelerating 
the convergence of certain continued fraction representations which are used 
in numerical algorithms; for instance, the algorithms in [lOllllj could be im- 
proved by using the bounds of sections 14.1.21 and 14.1.31 for accelerating the 
convergence. 



2 Qualitative behavior of Riccati equations 



In this section we deal with first order differential systems ([T]) with differen- 
tiable coefficients. We consider the ratio hn{x) = yn{x) /yn-i{x) satisfying 
the Riccati equation 



h'{x) = d{x) — {b{x) — a{x))h{x) — e{x)h{x)'^. 



(4) 



The label n, which is common for h and the coefficients a, b, d and e, has 
been dropped in @ for simplicity and because the analysis in this section 
is valid for any system, depending or not on a parameter n. The explicit 
dependence on n will be recovered in the next section. 
We have h'{x) = when h{x) = \^{x) with 




sign(e(x))i?(x) —r]{x) ± Y^ry(xp~+~s 
,ri{x) 



b{x) — a{x) 
2^\d{x)e{x)\' 



sign{d{x)e{x)), 



(5) 



We consider the case with real roots A^(x). Two distinct situations may 
occur: either d{x)e{x) > 0, or d{x)e{x) < but \r]{x)\ > 1. 

The condition d{x)e{x) > generally holds in the whole maximal interval 
of continuity of the functions because the coefficients d{x) and e{x) do not 
change sign under very general conditions (see, for instance, \12\ lemma 
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2.1]:^). Contrarily, when d{x)e{x) < the condition \rj{x)\ > 1 may hold only 
for a limited range of the variable x. In the first case {d{x)e{x) > 0) h{x) 
may have one zero or one singularity, but not both ( [12\ lemma 2.4]), while 
in the second h{x) may have both a zero and a singularity ( |13^ Theorem 
2.1]). We analyze the case d{x)e{x) > and assuming that no change of 
sign of h{x) occurs. For the case d{x)e{x) < 0, as the examples in section 
14.21 will show, similar arguments can be applied. 

In the sequel, we consider d{x)e{x) > 0. Without loss of generality, we 
take d{x) > 0, e(x) > and then A+(x) > and X^{x) < 0; if d{x) < 0, 
e{x) < we can consider the replacement y — )• —y or w ^ —w. In the next 
results, (a, h) is an interval where h{x) and the coefficients of the system are 
differentiable; a or 6 could be +oo or — oo. Depending on the value of h{x) at 
a"*" or different bounds can be established. First we consider h{a^) > 0. 
We enunciate three results and give a common proof. 

Lemma 1. If h{a^) > then h{x) > in {a,b) 

Theorem 1. Ifh{a+) > 0, A+(a;) is monotonic and h' {a'^)X^' {a^) > then 
{h{x) - \+{x))\+'{x) < in (a, 6). 

Theorem 2. // h{a^) > 0, A^(x) is monotonic and ft,'(a+)A^'(a+) < 
then either h{x) reaches one relative extremum at x^ S (a, b) (a minimum if 
A^'(x) > and a maximum if X^'(x) < 0) or {h{x) — A+(x))A+'(x) > in 
(a, 6). 

Proof. If /i(a^) > 0, then h{x) can not change sign continuously: it can not 
become zero because h'{x) > if < h{x) < X^{x). On the other hand, 
it can not change sign discontinuously; for this, starting with h{a'^) > 0, 
a value Xoo G {cL,b) should exist such that h{x^) = +oo but this is not 
possible because h'[x) < if h{x) > A+(x). 

Now, we consider that A"*"(x) is monotonic. We take the case A"'"'(x) > 0; 
the case X~^'{x) < is analogous. 

^AU that is required is that the system is satisfied by two independent sets of functions 



5 



Assume first that h'{a^) > 0; using ^ this means that < h{a~^) < 
A+(a+). And then, necessarily h{x) < A+(x) in (a, 6). Indeed, because 
A"'"(x) is monotonically increasing and the graph of h{x) is below the graph 
of X^{x) close to X = a, the graph of h{x) may touch the graph of A^(x) 
at X = Xg only if the first one has a larger slope at Xe, that is, if h'{xe) > 
A+'(xe) > 0; but if h{xe) = A+(xe) then h'{xe) = 0. 

If, contrarily, h'{a^) < then the graph of h{x) lies above the graph of 
A"^(x) close to X = a and there are two possibilities: either it remains above 
A"''(x) in all the interval or there is a point Xe G (a, b) where h{xe) = A'*"(xe) 
and h'{xe) = 0. Then the graph of h{x) crosses the graph of A^(x), which 
is an increasing function, and h'{x) > for all x > Xe- Therefore there is a 
minimum at Xg. □ 

Figure 1 illustrates the situations described in Theorems [1] and [2l 



T ' r 
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Figure 1: The characteristic root A"'"(a;) divides the plane in two regions: h'{x) > if 
< h{x) < A^(a;) and h'{x) < if h{x) > \^{x). The graph of hi{x) corresponds to 
the situation described in Theorem [1] while h2(x) corresponds to Theorem [2] when an 
extremum is reached. 

If, differently from theorems [T] and EJ we have h{a'^) < then h{x) may 
change sign once. But if it does not change sign and h{b^) < we are in 
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the previous situation. Indeed, with the change of variable x — )• —x and the 
change of function w{x) — )• —w{x), we have that the new ratio of functions 
h{x) = —y{—x)/w{—x) is such that h{a'^) > and the previous results hold 
in the interval [a,/3] = [—b,—a]. Then, we can write a common result for 
both cases. We only give the result corresponding to Theorem [TJ 

Theorem 3. Let h{x) he a solution of ^ with continuous coefficients and 
d{x) > 0, e{x) > 0. Suppose that either h{a^) > or that h{b^) < and 
take s = +, c = a+ in the first case and s = —, c = in the second. Then, 
h{x) does not change sign in {a,b), and if the characteristic root A'^(x) is 
monotonic and X'"{c)h'{c) > then 

i\h{x)\-\X'{x)\)^<Oyxe{a,b) 
ax 

Remark 1. The condition X'^'{c)h'{c) > is equivalent to 

(IMc)l-IA^(c)l)— (c)<0 

3 Bounds for first order DDEs 

Now, consider a first order difference-differential equation ([T]) and assume 
it holds for n > uq and that that the shift n — )• n + 1 is possible (true for 
continuous dependence on the parameter n). Then the solutions of ([T]) are 
also solutions of a three-term recurrence relation 

en+iy„+i(x) + (bn+iix) - an{x))yn{x) - dnVn-iix) = 0. (6) 

As in the previous section, we assume dn{x)en{x) > 0. 
Let be the roots of the algebraic equation 

Cn+lA^ + {bn+1 - an)Xn " = 0, (7) 

that is: 



A± = RnEni-Vn ± + 

Rn = \J dn/eni En = \J Cn/ En+l, ??n = ipn+l — «n)/(2 Y^^^nAi+l; 



(8) 
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If lim„_j.+oo / then lim„_i._|_oo | / | / 1, and if the coefficients 
are of algebraic growth as a function of n, Perron-Kreuser theorem (see [14^ 
Thm 4.5]) states that independent pairs of solutions {y^f}\y^^^} exist such 
that 

lim ^4t^ = 1, hm J-^ = l. (9) 

If f/„ > the minimal solution is yli^ and yi^^ is dominant, and therefore 
liuin^^oo Vn^ / Vn^ = . If ry„ < the roles are reversed. In both cases we 
have, for sufficiently large n, y„^^ y'n^ > and y^^l^yi?^ < 0. 

Remark 2. The minimal solution satisfies fjnyn/yn-i > for large n, while 
the dominant solutions are such that f]nyn/yn-i < for large n. 

Notice that the roots ([8]) are closely related to the characteristic roots 
of the Riccati equation ([5]): 



A^W = j^A-Vn{x) ± Vl + Vn{xn Vn{x) = ""^f^ - (10) 

V ^'^y^> 2^^ dn[x)en{x) 

As we have shown in the previous section, when A^(x) are monotonic they 
provide bounds for some solutions. On the other hand, if lim„^+oo A^/A^ = 
1 the function ratios have these bounds as limits. This explains why the 
bounds (fTOj) tend to be sharper as n becomes larger. Because of this, we 
refer to these bounds as Perron-Kreuser bounds. 

In section 13.21 we will obtain additional upper and lower sharp bounds 
starting from the bounds of Theorem [1] and using the three-term recurrence. 

Before this, it is important to stress that for the Perron-Kreuser bounds 
to hold, it is crucial that the characteristic roots are monotonic as a function 
of X. This, however, is a quite general situation, as we next see. 

3.1 Monotonicity of the characteristic roots 

The next result relates the monotonicity properties of the characteristic 
roots with monotonicity properties as a function of n which are known to 
hold for a large set of functions. 
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Theorem 4. Let yk{x), k = n,n — 1, be solutions of second order ODEs 
yki^)^^k{x)y'f,{x) + Ak{x)yk{x) = 0, withAk{x), Bk{x) continuous in {a,b) 
and Bn{x) = Bn-i{x). Assume that yn{x) and yn-i{x) satisfy a system ^ 
with dn{x)en{x) > and differentiahle coefficients. Then, en{x)/dn{x) is 
constant as a function of x, and if An{x) ^ An-i{x) the characteristic roots 
AJ(x) i fTOj) are monotonic in {a,b). Furthermore, dX^{x)/dx has the same 
sign as An-i{x) — An{x) and — r/^(x). 

Proof. Differentiating the first equation of the system ([T]) and ehminating 
yn-i and proceeding similarly with the second equation we have 

y'^{x) + Bk{x)y'^{x) + Akix)ykix) = 0, k = n,n-l, (11) 
with coefficients satisfying: 



BJx) - Bn-iix) = ^^44 - ^f¥\- 



(12) 



Now, because we are assuming that Bn{x) = Bn-i{x) the first equation im- 
plies that dn{x)/en{x) does not depend on x. Therefore, from the expression 
of the characteristic roots ([8]) we see that dX^{x)/dx has the same sign as 
—r]'^{x). All that remains to be proved is that An{x) — An-i{x) has the 
same sign as r]'^{x). But considering the second equation of (fT2]) and using 
that d'n{x)/dn{x) = e'j^{x) / en{x) one readily sees that An{x) — An-i{x) = 
'^\/dnix)enix)r]'^{x), which proves the theorem. 

□ 

Remark 3. If en{x)dn{x) < and rjnix)"^ > 1, it is also true that both roots 
are monotonic if Bn{x) = Bn~i{x) and An{x) / An-i{x), but X^{x)X^{x) > 
and X^'{x)X~'{x) < in this case. 

The case described in Theorem [J] is, for instance, the situation for Bessel 
functions, parabolic cylinder functions and the classical orthogonal polyno- 
mials when n is the degree of the polynomials. 
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3.2 Perron-Kreuser bounds 

In the following, we assume that rjn{x), fjnix), dn{x), e„(x) and hn{x) = 
yn{x) / yn-i{x) do not change sign for large enough n (say n > no). Notice 
that the sign condition for /i„ (x) is satisfied for large enough n when Perron- 
Kreuser theorem holds. An immediate application of Theorem [3] gives: 

Theorem 5 (First Perron-Kreuser bound). Let dn{x) > 0, e„(x) > and 
hn{x) = yn{x) / yn-i{x) with constant sign for n > no and for any x G ia,b). 
Let s = sign{hn{x)) and A^(x) as in Eq. [Td\) . Then, if hn{a'^) > and 
/i^(a+)A^'(a+) > or hn{b^) > and /i^(6")A^'(6") > the following holds 
in (a, b): 

{\hn{x)\ - F^{x)) A^'(x) < 0, n > no (13) 
F^{x) = Rn{x)(-srj^{x) + Vl + Vnixf) = , ^ ^ , (14) 

Further bounds can be obtained by iteration of ([6]), which we write: 

- d„ 6„+l - On + en+l ;r— (15) 



yn-l[x) V yn{.x) 

it is clear that for minimal solutions {fin{x)yn{x) / yn~i{x) > for large n), 
by substituting y„_|_i(x)/y„(x) by a lower (upper) bound we get an upper 
(lower) bound for 2/„(x)/y„_i(x). We only give the first iteration. 

Theorem 6 (Second Perron-Kreuser bound for minimal solutions). Under 
the conditions of Theorem\^ and if sfin > 0, s = sign{hn) then 

(|/i„(x)|-5^)A^'(x)>0, n>no (16) 



where 



gs+ ^ DARn ^ ^^^^ 

s{2Dnfin - r]n+l) + + ??n+l 



Dn = y/dn/dn+i; En, Rn and rjn given by (0) and ry„ by ^^). 
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The second superscript of the notation stands for the sign oifjnyn/Un-i- 

Notice that theorem [6] may be true for n = uq — I too, because Theorem 
Ois used in the proof with the shift n — )■ n + 1. 

The similarity of the second expression of ()14p with (jl7p indicates that for 
coefficients of algebraic growth we will generally have lim^^+co Fnl Sti' — ^• 

Further iterations are possible and this gives a convergent sequence of 
upper and lower bounds under the conditions of Theorem [5] and [6] and pro- 
vided that Perron-Kreuser theorem holds (which implies that the recurrence 
admits a minimal solutions). We don't prove this result, but the convergence 
of the sequence of bounds for the minimal solution follows immediately by 
using the same arguments considered in [Ij for the case of Modified Bessel 
functions of the first kind. 

We can also obtain additional bounds for dominant solutions by writing 

yn{z) 6„-a„_i dn-iyn-2{z) 

T\ ~ ^ TT ^^°> 

yn-l{x) Cn en yn-l{z) 

Differently from the case of minimal solutions, the sequence of bounds is not 
a convergent sequence. We give an explicit formula for the first iteration: 

Theorem 7 (Second Perron-Kreuser bound for dominant solutions). Under 
the conditions of Theorem\B\ and if sfin-i < 0, s = sign{hn), 

{\hn{x)\ - 5r)A^'(x) > 0, n > no + 1 (19) 

where 

Notice that the previous theorem can only be guaranteed to be true for 
n = no + 1, because Theorem [5] is used in the proof with the shift n — )• n — 1. 

The similarity of the first expression in ()14p with (j20p is clear. For 
coefficients of algebraic growth we will generally have lim„^_|_oo F^/S^^ = 1. 
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3.3 Turan-type inequalities 

Because upper and lower bounds are available for both when y„ 

is a minimal or a dominant solution (Theorems [5l [6] and [T]), upper and 
lower bounds for \yn/ yn-iWUn/ yn+i\ become available. The modulus can be 
skipped if yn/Vn-i does not change sign (as assumed earlier). With this: 

In < Ln{x) < ^"^f , < Unix) < Un, (21) 

yn+i{x) yn-i{x) 

where In = mina;{L„(x)} and Un = max2.{C/„(x)}. Many new Turan-type 
inequalities are found in section [J] by using this simple idea. 

3.4 Bounds of Liouville-Green type 

Using the difference-differential system ([T]) and the Perron-Kreuser bounds, 
bounds on the logarithmic derivatives can be established. We give the 
bounds obtained from the first Perron-Kreuser bound. 

Theorem 8. Under the hypothesis of Theorem\^ and if d\^/dx > Q (s = 
sign{yn{x)/yn~i{x))): 

s < s h V«n(2;)en(a;) V 1 + ^n\xY < s — -- (22) 

yn~l{x) 2 yn{x) 

IfdXf^/dx < the inequalities are reversed 

Two consequences follow. First, we observe that the ratios y'j^{x)/yk{x) 
are monotonic as a function of the discrete variable k. Second, because we 
are assuming that the shift n — )■ n + 1 is possible, we have both an upper and 
a lower bound for y'n/Vn- Upper and lower bounds could also be obtained 
by considering both the first and second Perron-Kreuser bounds. 

In the examples we will see that these bounds, after integrating the 
logarithmic derivative, are related to the Liouville-Green approximation 
for solutions of second order ODEs. In fact, using this analysis and by 
Liouville-transforming the first order system associated to the ODE y"{x) + 
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A{x)y{x) = 0, conditions can be established under which the LG approxi- 
mation for the solutions the ODE y" (x) + A{x)y{x) = are bounds for some 
of the solutions. We leave this analysis for a future paper. 

4 Applications 

We give a number of examples of application of the techniques described 
in the paper. We concentrate mainly on the case d„(x)e„(x) > 0. We also 
give two examples of application for monotonic systems with (i„(x)e„(x) < 0. 
The examples given by no means exhaust the functions for which the analysis 
is possible. 

4.1 Cases with dn{x)en{x) > 

We give examples which include classical orthogonal polynomials outside 
their interval of orthogonality. In all cases except the last one, Theorem H] 
holds. The last case is that of Laguerre functions of negative argument, for 
which Theorem [J] can not be applied but the characteristic roots are still 
monotonic and the same analysis is therefore possible. Some monotonicity 
properties for the determinants of some of these functions (modified Bessel 
functions, Hermite polynomials of imaginary order and Laguerre polynomi- 
als of negative argument) were considered in |15] . 

4.1.1 Modified Bessel functions 

These are solutions of x'^y" + xy' — (x^ + = 0. This was the case 
considered in detail in [1], and most of the results obtained in that paper 
are direct consequences of the more general results of the present one. 

4.1.2 Parabolic cylinder functions 

The parabolic cylinder function U{n, x) is a solution of the differential equa- 
tion y"{x) — + n)y{x) = 0, with coefficient A{x) = — + n) de- 
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pending monotonically on the parameter n (Theorem [J] holds). 

Considering the DDE satisfied by U{n,x) \16\ 12.8.2-3] and defining 

Unix) = e™'^U{n, x) we have: 



y'n-li^) = -^yn-l{x) + (n - l/2)yn{x). 

where n will be real and positive. For this system 



(23) 



"^^^^^ = ~ „ / ^n{x) = 77n+l(x), A^(x) = , (24) 

2Vra-l/2 rc=FV4n-2 + x2 

From [16, 12.9.1] we have /i„(+oo) = 0~ and /i[j(+cxd) = 0+ and because 
A~(+oo) = 0^ then theorem [3] holds, as well as theorems [5] and [6l Therefore 

Theorem 9. For n > 1/2 and x > the following holds 

2 U(n,x) 2 , , 

< ., , < , (25) 



X + y/An + 2 + x2 U{n-l,x) x + \/4n - 2 + x" 
The lower bound also holds ifn£ ( — 1/2, 1/2) and it turns to an equality 
ifn = -1/2. 

The lower bound is obtained from the upper bound and the application 
of the three-term recurrence relation: if Bfn(n, x) is a positive upper (lower) 
bound for U{n, x)/U{n — l,x), x > 0, then 

B^+i{n,x) = l/(x + (n + l/2)B„(n + l,x)) (26) 

is a lower (upper) bound for the same ratio. The process can be continued 
as m — )• +00 and the sequence is convergent (because U{n, x) is minimal). 

Now, consider ynix) = U{n,—x), which is also solution of ([23]) . Using 
the values of C/(n,0) and U'{n,0) [HI 12.2.6-7] it is easy to prove that 
/i„(0+) > 0, h'^{0+) > 0, n > 1/2, X > and then h'^{0+)dX+{0+)/dx > 
and Theorem [1] holds. The corresponding Perron-Kreuser bounds (theorems 
[5]and[7|), give: 



•^It is not important that the new functions are complex, because we are deahng with 
ratios; an alternative definition could be yn{x) = {—l)'^"^U(n,x). 
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Theorem 10. For n > 3/2 and x > the following holds 

X + \/4n — 6 + U{n, —x) x + \/4n — 2 + 

2n - 1 ^ C/(n- 1,-x) ^ 2n - 1 ^ ^ 

T/ie upper bound is also valid if n G (1/2,3/2). 



The upper bound in (j25|) has the same expression as (|27|l but with x in 
replaced by — x. Therefore: 

Remark 4. Theorems\^ and\T^ hold for all real x, hut for x < the lower 
hound of Theorem{^ only holds for all x < if n > 1/2. The lower hounds 
are sharper when x > 0. 

The fohowing Turan-type inequalities follows from Theorems [9] and [TOl 

Theorem 11. Let F{x) = U{n, xf/{U{n - 1, x)U{n + 1, x)). 
The following holds for all real x: 



n-3/2 n-1/2 , , , , n + 3/2 , , 

The first inequality holds for n > 3/2 and the rest for n > 1/2. For x < 
the third inequality also holds if n £ (—1/2,1/2). 

Finally, considering Theorem[8]and writing together the results for U{n, x) 
and U{n, — x) we have the next result. 

Theorem 12. For all real x and n > 1/2 the following holds: 

U'{n, x) 



- Vx74 + n + l/2< " ; '7 < -VxV4 + n-l/2 (29) 

U [n, xj 

The left inequality also holds for n > —1/2. 

These type of bounds are useful for studying the attainable accuracy of 
methods for computing the functions. In [17J, the following estimation for 
large x and/or n was considered for the condition number with respect to 
x: 

C^(C/(a,x)) = |xC/'(a,x)/C/(a,x)| ~ x^Jx'^/A + a, (30) 
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and similarly for V{a, x). The bounds (j29p prove that this a good estimation 
because it lies between the upper and lower bounds. From the previous 
discussion on the V{a, x) function, one can prove that similar bounds are 
valid for moderate x (x > 1 is enough); we consider later this function. 
Integrating (p9|) we have 

i^n+l/2(^)/i^n+l/2(y)<i3j^<i^n-l/2(x)/F„-l/2(^^ , ^ 

/ ^K 'f) „ (31) 

Fa{x) = exp UyGF/iTa) (x + 2^x^/4 + aj 

and, in particular, 

^ , . U(a,x) ^ , , , , 

F,+i/2(x) < < F„_i/2(x) (32) 

where 

F„(x)=expl-^W:-+a| I __ + W I- + 1 I (33) 




The bounds (j32p are useful for obtaining the range of parameters for 
which function values are computable within the arithmetic capabilities of 
a computer (overflow and underflow limits). These results confirms the 
estimations based on the Liouville-Green approximation used in |18j . 



o Iterated coerror functions and Mill's ratio: In particular, consider- 
ing Theorem[9]and the relation of parabolic cylinder functions U{n + 1/2, x) 
with the iterated coerror functions z"'erfc(x) [16t 12.7.7], n G N, the following 
follows: 

M„+i(x) < ^_^^^';(^) < M„(x), n= l,2,...;M„(x) = {x + ^ 2n + x"^)-^ . 
I erfc(xj 

(34) 



These inequalities appear in [19] . 

Theorem [9] also gives bounds on Mill's ratio {n = 1/2). From lower 
bound in Theorem ([9]) and the upper bound obtained by iterating with (p6|) 
we have 
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Theorem 13. Let r{x) = e^^/^ J^°° e ^^^^dt, then 



2 4 

< r{x) < (35) 



The lower bound was obtained in [20] and the upper bound in ^J. In 
our case, these results follow from a more general result. See also [22] for an 
alternative proof. 

Further iterations (see (p6]) ) give additional sharper bounds: 

Theorem 14. 

R2k+i <r{x)< R2k{x) (36) 
112 n 



Rn{x) = — — — ■ ■ ■ TfTT^, Tn = {x+ V^n + x^)/2 (37) 
X+ X+ X+ Tn{x) 

where, as usual we denote . . . = l/(a + 1/(6 + . . .)) 

o Hermite polynomials of imaginary variable A similar analysis to 
that for U{n,—x) can be carried for the PCF V{n,x). Indeed, yn{x) = 
y(n,x)/r(n + 1/2) is a solution of ([23|) and hn{x) = yn{x) /yn-i{x) is such 
that /i„(0^) > 0. Two situations take place depending on the values of n. 
First, if n G {2k - 1,2/c), /c G N, then /i^(0+) > and the upper bound of 
Theorem [10] holds in this case and for all x > while the lower bound will 
hold for n G (2A:, 2A; + 1). Contrarily, if n G (2A;, 2A; + 1) then h'^{'d+) < 0, 
while hn{+oo) > 0, and the upper bound only holds for large enough x; a 
similar situations occurs with the lower bound when n G {2k — 1,2k). 

We only consider the first case. Then, using the relation of V{n-\-l/2, x), 
n G N, with Hermite polynomials [16, 12.7.3] we get: 

Theorem 15. 



< " + ^'y + "' ,x>0,nG(2fe-l,2fc),fcGN (38) 



V{n-l,x) 2 



X + v 4n — 6 + V{n,x) . , , x , , s 

X > 0, nG 2fc,2fc + l , fc G N 39 

2 K (n — 1, xj 
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. H2k+i{ix) 



<x + \J^k + 2 + x2, X > 0, A; = 0, 1, 2 . . . (40) 



. H2k-i{tx) ^ ^ ^4fe _ 2 ^ 3,2)-!^ X > 0, fc e N (41) 



H,U.^ >,tlV|,,,M,..M. (42) 



H2k-i{ix)H2k+i{'ix) y A; + 1/2 

Hermite polynomials of imaginary argument were also considered in 
|15j . The well-known Tm'an-type inequality for Hermite polynomials [23] 
Hn{x)'^ — Hn-i{x)Hn+i{x) > 0, X G M, does not hold on the imaginary axis, 
but a similar property Hn{ix)'^ — ^y {n — l)/{n + l)Hn-i{ix)Hn+i{ix) > 
holds true for all x > if n is even. 

4.1.3 Oblate Legendre functions 

These are Legendre functions of imaginary argument, which are functions 
appearing in the solution of Dirichlet problems in oblate spheroidal coordi- 
nates [TU]. Denoting 

and using the differential relations [161 14.10.4-5] we have 

P'n{x) = ^ 2 {nxpu{x) + (n + m)pi,_i(x)} 

1 + X (44) 

p'n-i{x) = , 2 {-nxpj,_i(x) + (n - m)p^(x)} 

1 + X 

and qn{x) = Q^{ix), Q"^ being the second kind Legendre function, satisfies 
the same system. We consider n > m and x > 0. This is again an example 
for which Theorem H] holds. The roles played in this case by the functions 
Q^{ix) and P^(ix) are very similar to the roles of U{n,x) and V{n,x) in 
the previous section. We omit details and only summarize the main results. 

Theorem 16. The following holds for x > and real n > m > 



^ i(fx) ^ n 




-1 



n + m 
< \ 45 
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Q'!^_i{ix) 



> 



n + m 



(46) 



nx + (n + 1)\ 1 + x"^ 



m 



1 < 



n + m + 1 



< 



'(n + 2)2 -m^ 

2 2 

n — m 



Theorem 17. T/ie following holds for x > and integer n,m, n > m: 



(47) 



< -I " \. \ < 



n 



P^_i{ix) n-m 



m 



X + \ 1 + x^ ^ 



n 



, n — m odd (48) 



n — m 



nx + {n — l)i 1 + x'^ 



m 



(n-1) 



< 



p 



n — m even 



n—l 



[IX 



p:rizx) 



< 1 + 



1 



n — m 



-, n — m odd 



(49) 
(50) 



p:r_,iix)p:r+iiix) 

For m = we have Legendre polynomials. If n is odd, we have Pn{ix)'^ < 
and therefore Pn{ixY — (1 + ^/T^)Pn-i{ix)Pn+i{ix) > 0. It appears, as 
numerical experiments show, that in this case the same Turan inequality 
that holds in the real interval (—1,1) [24J also holds in the imaginary axis 
if n is odd: Pn{ixY — Pn-iiix)Pn+iiix) > 0; the same is not true if m / 0. 



4.1.4 Laguerre functions of negative argument 

Next we consider an example for which Theorem S] can not be applied but 
the analysis is possible because the characteristic roots are monotonic. 

Consider the Laguerre functions yu,a{x) = L"(— x), x > 0. Using well 
known recurrences and differentiation formulas, we have 

yl+l,a-li^) = yv,a{x) ^^^^ 

Xy'p,a{x) = + x)yi,^a{x) + {y + l)?/i.+l,a_l 

and 

{v + l)yy+i^a-i{x) = (a + x)yu,aix) + xyu-i,a+i (52) 
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Considering \25\ Theorem 2] it follows that y^^a is a dominant solution of 
the recurrence (|52p in the direction of increasing (and decreasing a). 

With h{x) = yu,a{x) / yu+i,a-i{x) , the positive characteristic root A+(x) 
of the associated Riccati equation turns out to be increasing if u > —1 and 
a > 0. On the other hand, it is easy to check that for these values /i(0"*") > 
and /i'(O^) > 0. Theorem [1] holds and A^(x) is a bound: 

Theorem 18. For any a > 0, v > —1 and x > the following holds 



X 



On the other hand, from the recurrence <\52h we have 
L^(-x) fa + x , X Klli-x)^ 



+ : (54) 



and from this we obtain the second Perron-Kreuser bound: 
Theorem 19. For any a > —1, v > Q and x > the following holds 
^°+i^(-x) g + X - 1 + V(a + X+IY + 4.UX 

And from these bounds we get the following Turan-type inequalities: 

Theorem 20. For any > and a > 0, x > the following holds: 
V a ^ LZ-\^-x)LZt\{-x) ^ y 



u + la + 1 L'^i-x) L'^i-x) u + l 
A second independent solution of (j5ip which is a minimal solution of 

(|52p as z^ — )■ +00 follows from \2h\ Theorem 2]. Bounds can be also obtained 

for this solution. We omit the details. 

Other bounds and inequalities can be obtained using other recursions 

or using relations between contiguous functions. For example, using 

18.9.13], we have: 

LZ{x) ^ LZ{x) ^ ^ 

and upper and lower bounds for L°(— x)/L°_^(— x) follow from the previous 
results. As a consequence of this new bounds, one can prove the following 
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Theorem 21. 



u+1 



V 



< 



< 



u+lu+a— 1 



(58) 



where the first inequality holds for u > a > — 1 and the second for u > 
u + a > 1. 

For positive x, it is known that L"_^(x)L"_^]^(x)/L"(x)^ < 1 [23]. For 
negative argument we have an upper bound greater that 1, which suggests 
that the Turan-type inequaUty for positive x does not hold for negative x, 
as numerical experiments show. 

4.2 Two examples with (i„(x)e„(x) < 

The DDEs corresponding to a pair {p„(x),p„„i(x)} of classical orthogonal 
polynomials satisfy (i„(x)e„(x) < in their interval of orthogonality because 
this is a necessary condition for oscillation \\.2\ Lemma 2.4]. However, for 
values of the variable for which the polynomials are free of zeros, one can 
expect that r]n{x)'^ > 1 and that the DDE becomes monotonic (r/„(x)^ < 1 
is also a necessary condition for oscillation [13, Thm. 2.1]). This is the case 
of Laguerre and Hermite polynomials for large enough x > 0. We consider 
these two examples. 

4.2.1 Hermite polynomials 

Hermite polynomials satisfy 



H^{x) = 2nF„_i(x), 
H'^_^{x) = 2xHn-i - Hn{x) 



(59) 



We have r]n{x) = x/\/2n and r]n{x) > 1 if x > \/2n (monotonic case). The 
characteristic roots are both of them positive 



AJ(x) = X ± \J x'^ — 2n. 



(60) 
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Defining hn{x) = Hn{x) / Hn^i{x) we have that /i„(+oo) = +00 and 
h',^[+co) > because the coefficient of degree n of Hn{x) is positive. Then 
hnix) > Xni^) enough x > because h'^{x) > only if hn{x) < A~(x) 
or hn{x) > A^(x), but /i„(+oo) > A~(+oo) = 0^. Then, we have that 
hn{x) > A+(x) for large x. And because A+'(x) > if x > \/2n, then, 
necessarily: 

hn(x)= , " A > X + -2n,x> V2^. (61) 
Hn-i{x) 

We can iterate the recurrence relation. Contrary to the case e„(x)(i„(x) > 0, 
we will not obtain sequences of lower and upper bounds, but only lower 
bounds. Writing 

hn+i{x) = 2x — 2n/hn{x) (62) 

and using (j6ip we get a lower bound for /i„+i(x). We shift the parameter n 
and get 

hn{x) > x+ ^/x^ - 2{n - 1), X > V2(n- 1). (63) 

This improves Eq. (j6ip and enlarges the range of validity of the bound with 
respect to x, but reduces the range of validity with respect to n (n > 2). 
The next iteration gives a bound for n > 3: 

hn{x)> F{n,x),x > ^2(n-2) 

F{n,x) = (n - 2)-i[(n - 3)x + (n - l)v^x2 - 2(n - 2)] 

Taking into account the largest zero of Hn{x) is larger than the largest 
zero of Hn-i{x) the previous bounds give bounds on the largest zero of 
Hn{x). We see that the largest zero of Hn{x) is smaller than \j2[n — k) if 
n > k. 

We consider just one more iteration and get 

hnix) > 2x-2(n- 1)/F(n- l,x) = G{n,x), x > ^J2{n - 3) (65) 

and if G(n, ^2(n - 3)) > then G(n, x) > if x > ^/%n^^, and the 
largest zero will be smaller than \j2{n — 3); this condition is met if n > 7. 
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A sharper bound has recently appeared in the Uterature |26j vahd for all n. 
However, the result is sharper than previous results, like for instance those 
in [27], which is interesting given the simplicity of the analysis. This reflects 
the fact that the bounds on function ratios (our main topic) are sharp. 

4.2.2 Laguerre polynomials 

We give some results for Laguerre polynomials omitting details. Defining 
= — L"(x)/L"_]^(x), we have /i"(+oo) = +00 and /i"'(+oo) = +00 
and, proceeding similarly as before: 



2nh^{x) > X - {2n + a) + y^{x - 2n - a)^ - ^n{n + a), 
X >2n + a + 2\Jn[n + a) 

and after the first iteration of the recurrence we have: 



(66) 



2n/i^(x) > /(x), X > 2n* + a + 2^n*(n* +a), n* = n-l, 



(67) 



f{x) = X — {2n + a) + y^(x — 2n* — a)^ — 4n* (n* + a) . 

This proves that the largest zero of L"{x) is smaller than x* = 2n + 
a — 2 + Y^(n — l)(n — 1 + a), provided that f{x*) > 0, which is true if 
a > {n — 1)^^ — (n — 1), n > 2; notice that values a < —1 are allowed for 
large enough n. The bound in [28] is slightly sharper, and is improved in 
[26]. 

Further iterations are possible, but not so easy to analyze. The next 
iteration will give a bound 

2nh^{x) > g{x), x > 2{n - 2) + a + 2y/ {n - 2){n - 2 + a) = x* (68) 

X* is an upper bound for the largest zero provided that g{x*) > 0. This 
condition is met for a larger a range as n becomes larger. For n > 10, this 
holds for any a > —1. The bound (j68p is of more limited in terms of n but 
numerical experiments show that it is sharper than the bound in [26] for 
a < 12 
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We expect that lower bounds for the smahest zero can be also obtained 
with a similar analysis. 

The main message, as before, is that the bounds on function ratios are 
sharp for large x because they give the correct asymptotic behavior as x ^ 
+00, but also for moderate x given the sharpness on the bounds on the 
largest zero. 
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